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Abstract 

The physical cluster originally conceived by Hill (J. Chem. Phys, 23, 617) is generalized for 
the case > 2 in a novel way. Contrary to Hill's pairwise generalization, this definition assures 
that all constituent molecules of the cluster have insufficient kinetic energy to escape, and avoids 
the spurious implication of Hill's generalization that the molecular velocities are to be increasingly 
equal. 
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I. INTRODUCTION 



The rate of a phase transition obviously depends to a large extent on the rate at which the 
molecules can rearrange themselves, and in the case of condensation especially, this rate is 
much smaller than the typical rate of molecular translation. A condensation process necessi- 
tates diffusion along the concentration gradient and would be thermodynamically forbidden 
were it not for the decrease of chemical potential associated with the transformation, due to 
the attractive interactions dominating in the condensed phase. However, in the very early 
stage of the phase transformation, insufficient condensed phase is present to afford this sta- 
bilization. In these cases, the transport of matter must occur along the chemical potential 
gradient. The great reluctancy by which this fluctuation takes place is one bottleneck in 
the condensation process and was first identified by Gibbs, who expressed the free energy 
of a small cluster as the sum of a surface and a bulk term, taking the surface term to be 
proportional to the specific surface free energy of the bulk liquid This approach^was 
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further extended into a theory for the kinetics of the process by a number of authors 
This classical nucleation theory, as it has since become known, predicts a nucleation rate in 
fortuitous accord with experiment for water condensation in a narrow temperature interval. 

In efforts to side track the assumptions of bulk properties for small clusters on the order 
of 50-100 monomers, much computer simulation, following the pioneering work of Lee et 
al. [5|], has been carried out over the years, and new theories developed. However, the 
concept of a cluster has in these simulations often been quite arbitrary. Most often based 
on a simple distance criterion, as in e. g. the Stillinger \&\ or LBA 5[ cluster, there is no 
compelling argument a priori for the particular dimension of the cluster. The same problem 
of arbitrariness is faced when exchanging the distance criterion for one of potential energy 

Q. 



Of particular elegance is then the definition due to Hill js], who regarded molecules in 
regions of phase space of negative total energy to form a cluster, in analogy with how "bound 
states" are usually defined in the rest of physics, for instance between atoms in molecules, 
or between celestial bodies and satellites in astronomy. Moreover, the concept of the critical 
temperature follows directly from this notion, as the effective temperature where the kinetic 
energy is so great that no bound states can be formed, even for the infinite cluster, i. e. bulk. 
Following Hill, we denote such a cluster, a "physical cluster." However, Hill's generalization 



of his cluster definition past tlie dimer is based on a pairwise expansion in which monomers 
are not bound to the cluster as a whole, but to each other, and while rigorous, is not 
completely intuitive. In fact, for large clusters, his definition entails that the velocities of 



the monomers become increasingly harmonized jo]. In the author's view, nevertheless, one 
of the most important concepts of the Hill cluster is its inherent non-locality. Such non- 
locality has been exploited in other theories of nucleation, the most prominent being the 
density-functional approach of Oxtoby and Evans [10]. 

Consider the set 3^ of N particles ordered from 1 to A^. The particles have positions 
{ri}i , velocities {vi}^ and masses {nii}^ . Let U{{ri}) denote the potential energy of 
the system and let U{{fi}) have the property that U{{fi}) — )■ 0, when the smallest mutual 
distance in {fi} approaches positive infinity. Let us first recall Hill's definition of the physical 
dimer. 

Definition 1. The two particles i and j constitute a physical dimer if \ii{vi — VjY < 
—U{fi,rj) where fi 



rriimj 
rrii +mj 



Let us now define two ordered subsets of D^v denoted by A and B such that AUB = Dat. 
We shall denote the mutual energy of interaction between these two subsets as U (A, B) for 
conciseness. Furthermore, we define to be the sum of the elements of {rrii} corresponding 
to the particles of A. Likewise, let denote the mass-weighted average velocity of the 
elements of set A. Our generalized definition of Hill's dimer to the cluster is given below. 

Definition 2. The ordered set Djv of N particles constitutes a physical cluster if for all sets 
A ^ ant/ B = ©jv \ A ^ 0, have that 1^am{va - v^f < -f/(A, B) where fim = 

II. PARTITION FUNCTION OF THE HILL ENSEMBLE 

To treat the physical cluster in the theory of statistical mechanics, we need to derive 
its partition function. For the sake of notational simplicity, we shall concern ourselves only 
with the cluster at rest {vcm = 0). Center-of-mass motion will be assumed separable from 
the internal degrees of freedom. Hill js] gives the partition function for the dimer and the 



explicit expression for the probability 



p[-(3Uif,,r,)] = erf ( ^-/3f/(rl, f,) ) - l=e^^^^^-^^^) ^ -(3U{f,,fj) (1) 



that two molecules a physical cluster at inverse temperature P = 1/kT, where k is Boltz- 
mann's constant and T is the absolute temperature. Because of the strong analogy between 
Definitions [T] and |21 we are motivated to apply this probability function to each term in 
Definition HI but care must be exerted so as to avoid double counting. We first consider 
the sets Aj that contain the only element molecule j = 1,2, . . . , N. The probability that 
molecule j is part of the cluster is 

p,=p[-f3U{A„3N\A,)] (2) 

Next we consider the sets Aj^ that contain only the elements j and k. The probability that 
the dimer j, k is part of the cluster is 

Pjk = p[-pU{Ajk, \ Ajk)] (3) 

We continue to define probabilities like these, over greater and greater subsets Ajk... until 
we reach the subset that has N/2 elements where we must stop to avoid double-counting. 

The partition function for the cluster now follows upon substituting the canonical prob- 
ability function e~^'^^^^^^'^ by e-~^^^^'''^^ Y[jPjY[j<kPjk ■ ■ ■ the definition of the partition 
function. In the quantum case, we must deal with the subtle issues raised by the Heisenberg 
uncertainty principle, but in the classical case, this immediately gives 

Qcm = / dlrlje-'^^"-" Up^Up^^ n P^-^' • • • (4) 

j j<k j<k<l 

This partition function in essence defines a new ensemble: the ensemble of the physical 
cluster, or "bound" states. Its probability distribution is not of Boltzmann form. We shall 
refer to this new ensemble as the "Hill ensemble" in his honor. 



III. CONCLUSION 



Definition [2] is very general and encompasses both clusters of molecules or atoms, or 
the molecules themselves. However, the physical cluster is not defined for all potentials. 
For instance, the harmonic spring does not satisfy the requirement that the potential energy 
vanish for large distances. In physical reality, however, such an unbounded potential is never 
encountered, and all molecules, and especially their clusters, are prone to disintegration at 
high enough temperature. 
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Because the products run over all dimers, trimers, . . ., A^/2-mers, the numerical com- 
plexity in verifying whether molecules constitute a cluster according to Definition [2] is 
(A^/2)!. Hence, numerical experiments on the physical cluster are prohibitively expensive 
beyond the very smallest ones. However, generally it is clear that for T — )■ 0, we have 
p[—(3U{-)] —7- 1. Therefore, the Hill ensemble will reduce to the canonical ensemble for 
low temperatures. This explains the observation by Lee, Barker and Abraham jsj that the 
thermodynamic properties of Lennard- Jones clusters are largely independent of the choice 
of cluster constraining radius at low temperature. 
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